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^ We show that in four or more spacetime dimensions, the Einstein equations for gravita- 

^ , tional perturbations of maximaUy symmetric vacuum black holes can be reduced to a single 

. second-order wave equation in a two-dimensional static spacetime, irrespective of the mode 

of perturbations. Our starting point is the gauge-invariant formalism for perturbations in 
an arbitrary number of dimensions developed by the present authors, and the variable for 
the final second-order master equation is given by a simple combination of gauge-invariant 
variables in this formalism. Our formulation applies to the case of non-vanishing as well as 
vanishing cosmological constant A. The sign of the sectional curvature K of each spatial 
section of equipotential surfaces is also kept general. In the four-dimensional Schwarzschild 
background with yl = and K = 1, the master equation for a scalar perturbation is identical 
to the Zerilli equation for the polar mode and the master equation for a vector perturba- 
tion is identical to the Regge- Wheeler equation for the axial mode. Furthermore, in the 
four-dimensional Schwarzschild-anti-de Sitter background with yl < and -ft' = 0, 1, our 
equation coincides with those recently derived by Cardoso and Lemos. As a simple appli- 
cation, we prove the perturbative stability and uniqueness of four-dimensional non-extremal 
spherically symmetric black holes for any A. We also point out that there exists no simple 
relation between scalar-type and vector-type perturbations in higher dimensions, unlike in 
four dimension. Although in the present paper we treat only the case in which the hori- 
^ ^ zon geometry is maximally symmetric, the final master equations are valid even when the 

^ \ horizon geometry is described by a generic Einstein manifold, if we employ an appropriate 

, reinterpretation of the curvature K and the eigenvalues for harmonic tensors. 
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§1. Introduction 

Perturbative analyses of four-dimensional black hole spacetimes have provided 
useful tools for the investigation of astrophysical problems, such as gravitational wave 
emission from gravitational collapse and black holesP''^''^ as well as for the investi- 
gation of more fundamental problems, such as the stability and uniqueness of black 
holes.CJ'EJ''^''^ In these analyses, a key role is played by the fact that the Einstein 
equations for perturbations can be reduced to a second-order ordinary differential 
equation (ODE) of the self-adjoint type through harmonic expansion. Such a master 
equation in the Schwarzschild background was first derived by Regge and Wheeler 
for the axial modejEJ^ and subsequently by Zerilli for the polar 

mode.injnni 

Later, 

the formulation was further extended to the Kerr background by Teukolsky.^^EU 
Recently, it was also extended to the Schwarzschild-de Sitter and Schwarzschild- 
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anti-de Sitter backgrounds by Cardoso and Lemo^^'^ and utilised to determine 
the frequencies and the behaviour of quasi-normal modes P^'^'^'^'^ 

Judging from this experience in the four-dimensional case, it is expected that an 
extension of the formalism to higher dimensions would be quite useful for the inves- 
tigation of fundamental problems as well as phenomenological problems in higher- 
dimensional gravity theories, which are currently subjects of intensive research. In 
particular, such an extension will provide an exact basic equation for the quasi- 
normal mode analysis of higher-dimensional black motivated by the 
AdS/CFT issueP 

Such extensions already exist in some limited cases. For example, the present 
authors developed a gauge-invariant formalism for perturbations in spacetimes of 
arbitrary dimension greater than three in a previous papeJ^ (which is referred to 
as KIS2000 in the present paper) by extending the formulation presented in Refs. 
I24() and I25|). There, assuming that the background spacetime of dimension n + m 
possesses a spatial isometry group G„ that is isomorphic to the isometry group 
of an n-dimensional space /C"" with constant sectional curvature K, we expanded 
perturbation variables in terms of harmonic tensors on /C" and expressed the per- 
turbed Einstein equations as a set of equations for gauge-invariant variables in the 
m-dimensional orbit space constructed from the expansion coefficients. In this for- 
malism, the gauge-invariant variables are grouped into three types, the tensor type, 
vector type, and scalar type, according to the type of harmonic tensors used to 
expand the perturbation variables. Each of these types of variables obeys an in- 
dependent closed set of equations. In particular, we showed that in the case that 
m = 2 and the background spacetime satisfies the vacuum Einstein equations, the 
perturbed Einstein equations for both a vector perturbation and a tensor perturba- 
tion reduce to a single second-order wave equation — which is called a master equation 
in the present paper — in the two-dimensional orbit spacetime spanned by the time 
and radial coordinates. *•* This result is independent of the sign of K and holds for 
both vanishing and non-vanishing cosmological constant A. 

For a four-dimensional Schwarzschild black hole background, a vector perturba- 
tion corresponds to the axial (odd) mode, and the master variable and the master 
equation presented in KIS2000 correspond directly to the Regge- Wheeler variable 
and the Regge- Wheeler equation, respectively. Hence, we already have an extension 
to higher dimensions of the Regge- Wheeler formalism for the axial mode. However, 
such an extension for a scalar perturbation, i.e., an extension to higher dimensions of 
the Zerilli formalism for the polar (even) mode in four dimensions, has not yet been 
obtained, although simple master equations have been derived for a scalar perturba- 
tion in constant curvature spacetimes using various methods .'^''^''^''^''^''^''^''^''^ 

The main purpose of the present paper is to construct such an extension of the 
Zerilli formalism to spacetimes of arbitrary dimension greater than three. To be 
precise, we show that under the assumption of the G„ symmetry of the background 

*•* The master equations for these types of perturbations were first derived by Mukohyama in 
Ref. I29II for the case of a maximally symmetric background. In KIS2000, it was shown that the 
same equations hold without the assumption of maximal symmetry if the background possesses G„ 
symmetry. 
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spacetime, the gauge-invariant vacuum Einstein equations given in KIS2000 for a 
scalar perturbation in a system of greater than three dimensions can be reduced to 
a single second-order wave equation in the two-dimensional orbit space for a master 
variable <P consisting of a simple combination of the gauge-invariant variables. This 
equation is identical to the Zerilli equation for the polar mode in the four-dimensional 
Schwarzschild background. By the generalised Birkhoff theorem, the symmetry as- 
sumption requires the background spacetime to be static, and its metric is uniquely 
determined by the normalized curvature K of symmetry orbits, the cosmological 
constant A, and a parameter M representing the black hole mass, provided that the 
spacetime is not of the Nariai type,!^ which is assumed in the present paper. In 
order to allow for the widest possible application of the present formulation, these 
parameters are kept arbitrary. For example, the formulation describes a scalar per- 
turbation of a higher-dimensional Schwarzschild black hole for K = 1 and yl = 0, 
while it describes a scalar perturbation of a Schwarzschild-anti-de Sitter spacetime 
for yl < and M 7^ 0. We also show that in the special case M = 0, the master 
equation derived in the present paper is essentially equivalent to the master equation 
derived by Mukohyama in Ref. I29|l . Further, for completeness, we show explicitly 
that vector and tensor perturbations obey master equations whose structures are 
the same as in the case of a scalar perturbation, and we give the corresponding 
effective potentials for these equations. We also point out that the simple relation 
between the polar mode and the axial mode, which was first obtained by Chan- 
drasekhar and DetweileJ^ for a four-dimensional Schwarzschild black hole and then 
extended to the four-dimensional asymptotically anti-de Sitter case by Cardoso and 
Lemos,^'^ does not hold in higher dimensions. 

The present paper is organized as follows. In the next section, we show that the 
gauge-invariant perturbation equations for a scalar perturbation given in KIS2000 
are equivalent to a set of first-order ODEs with a linear constraint for three gauge- 
invariant variables, after carrying out the Fourier transform with respect to the time 
coordinate. Then, in §3, this set of equations is reduced to a single master equation 
of the Zerilli type in the case of a non-vanishing frequency. We further show that this 
master equation can be converted into a wave equation for a master field ^ in the 
two-dimensional orbit space by the inverse Fourier transformation and that all gauge- 
invariant variables can be represented as combinations of this field and its derivatives. 
A subtle point regarding static perturbations in this derivation is treated in §4. In 
§5, we rewrite the master equations for tensor and vector perturbations given in 
KIS2000 in the same form as that for a scalar perturbation. Section 6 is devoted to 
discussion. There, using the master equations for four dimensions, the perturbative 
stability and uniqueness of the Schwarzschild-de Sitter and Schwarzschild-anti de- 
Sitter black holes in four dimensions is proved. We also discuss the relation between 
the scalar and vector master variables and the extension of the formulation to the 
case in which the horizon geometry is described by a generic Einstein metric. Because 
the calculations required to derive the master equations are quite lengthy, most of 
them were done by symbolic computation with Maple. 
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§2. Basic equations 

In this section, wc show that the gauge-invariant equations for scalar pertur- 
bations given in KIS2000 reduce to a set of first-order differential equations with a 
linear constraint for three gauge- invariant variables in a maximally symmetric black 
hole background. 

2.1. Perturbation equations 

The general form of the metric of a G„-symmetric background spacetime con- 
sidered in the present paper is given by 

ds'' = gab{y)dy''dy' + r\y)dal (2-1) 

where gab is the Lorentzian metric of the two-dimensional orbit spacetime, and da"^ = 
^ij{z)dz^dz^ is the metric of the n-dimensional G„-invariant base space /C" with 
normalized constant sectional curvature K = 0, ±1. Hence, the dimension of the 
whole spacetime is n -|- 2, and the spherically symmetric case corresponds to i^T = 1. 

Scalar perturbations of this spacetime can be expanded in terms of the harmonic 
functions S on /C", which satisfy the equation 

(A„ + fc2)s = 0, (2-2) 

where A„ is the Laplace-Beltrami operator on /C*, and is its eigenvalue. For 
K = 1, k^ takes discrete values, 

A;2 = /(/ + n-l), Z = 0,l,2,--- , (2-3) 

while for K < 0, A;^ can take any non-negative real value. Each harmonic mode of 
the metric perturbation SgMN can be written 

^gab = /afoS, Sgai = rfa^i, 6gij = 2r'^{HL^ijE> + Hx^ij), (2-4) 

where 

Si = -^AS, (2-5) 

^i3 = -^Dibjn+^-iijn, (2-6) 

and Di is the covariant derivative with respect to the metric on /C". Here and 
hereafter we omit the indices k labelling the eigenvalues of the harmonics and the 
summation over them. 

Note that the modes with fc = and k"^ = uK require special treatment. First, 
for the A; = mode, S is a constant, and and Sjj are not defined. This mode corre- 
sponds to a perturbation that possesses the spatial symmetry Gn of the background 
spacetime. Because the only perturbation with the spatial symmetry Gn satisfying 
the Einstein equations consists simply of a change of the background metric in the 
form of a shift of the mass parameter of the black hole, as a consequence of the 
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Birkhoff theorem, we do not have to consider this case. Next, k'^ = nK occurs only 
for K = 1, and this corresponds to / = 1. For such modes, vanishes, and the vari- 
able Ht is not defined. We show below that these modes have no physical degrees 
of freedom. 

For modes with k'^(k'^ — nK) ^ 0, the following combinations provide a basis for 
gauge-invariant variables for metric perturbations of the scalar typeP^ 

F = Hl + ^Ht + ^D\Xa, (2-7a) 

Fab = fab + DaXt + DbXa. (2-7b) 

Here, Da represents the covariant derivative with respect to the metric gab in the 
two-dimensional orbit space, and Xa is given by 

Xa = r (fa + rDaHT) • (2-8) 



k V-'" ' k 

Because the Einstein tensors Gmn have the same structure as the metric per- 
turbation (|2-4)) under the harmonic expansion, the vacuum Einstein equations for a 
scalar perturbation are given by the following set of equations: 

Eab = r^-^Eab = 0, Ea = r^-^Ea = 0, El= v^-^El = 0, Et = t'^-^Et = 0. 

(2-9) 

The explicit expressions of Eab, Ea, El and Et in terms of the gauge- invariant vari- 
ables (|2-7|) can be obtained by specializing the general formula given by (63)"(66) in 
KIS2000 to the case in which the orbit space is two-dimensional (m = 2). As these 
expressions are long, we give them in Appendix \K\ 

The equations in (|2-9j) are not independent, due to the Bianchi identities, which 
are written 

^Dair^E") - -El + Et = 0, 2-lOa 

r-^ r n kr 

\D,{r^&J + -Ea - u^El = 0. (2-lOb) 

Since we are only considering modes with /c^ > 0, it follows from these equations 
that El = and Dh{r^E^) = are automatically satisfied if Ea = Et = holds. 
From HA-lb|) and (|A-ld|) . the latter equations are equivalent to the following set of 
equations: 



F« = -2(n-2)F, (2-11) 
Dk{F',-2F6i) = 0, (2-12) 



where 

Fab = r^-^Fab, F = r'^-^F. (2-13) 

Note here that the equation Et = is not obtained for the mode Z = 1 in 
the spherically symmetric case, because in that case Sjj vanishes. However, we can 
assume that this equation holds in this case too, by regarding it as a gauge condition, 
as is shown in Appendix IbI 
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2.2. Fourier transformation 

Because solutions to the vacuum Einstein equations with the spatial symmetry 
Gn are always static as asserted by the Birkhoff theorem, our background metric can 
take the form 

ds'' = -f{r)dt^ + /(^O ^ ^^'^^^ 

with 

f{r)=K-^-\r\ (2-15) 

provided that the spacetime is not of the Nariai type, which is assumed in the present 
paper. Here, A is related to the cosmological constant A by 

\= (2-16) 

n(n + l) ^ ' 

When the metric ()2-14|) with H2-15() describes a black hole spacetime, the G„-invariance 
of the spacetime implies the maximal symmetry of the spatial section of the event 
horizon. For this reason, we refer to the black hole described by the metric (|2-14() 
as a maximally symmetric black hole.*\ For = 1, A = 0, (|2-14|) is the stan- 
dard higher-dimensional Schwarzschild metric — also referred to as the Tangherlini 
metricP^ — with mass given in terms of the parameter M 

where An = 2^^''+'^^ ^ / r[{n + l)/2] is the area of a unit n-sphere and G denotes the 
n + 2-dimensional Newton constant. 

Because the non-vanishing connection coefficients for this metric are given by 

f f f f 

■pt J -pr J J 7-ir J /Q 1 Q\ 

^^"2/' * ~ ~ 2/' v^'-*-"/ 
the equations Dh{r'^E^) = can be written 

r'^dt&t + drir'^El) = 0, (2- 19a) 

r^dtEl + drir^K) + "^El - = 0, (2.19b) 

in the (t, r) coordinates. From the latter equation, we find that El = Q holds if 
the equations Db{r'^&^) = El = E; = are satisfied. Hence, taking account of the 
argument involving the Bianchi identities in the previous subsection, we only have 
to consider the equations (OTH) . (OT^ . E^ = and E^ = 0. 

*^ It is known that, in higher dimensions, the metric 12-141 1 with 12-15I I can describe many 
diflerent black hole solutions to the vacuum Einstein equations, with the simple replacement of da^^ 
by the metric for any Einstein manifoldP^^^ Among such solutions, the G„-symmetric metric is 
the maximally symmetric solution. Our formulation also holds in the case that da^ is given by a 
generic Einstein metric, as mentioned in §6. 
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Next, we reduce these equations to first-order ODEs for three variables. First, 
let us introduce the variables X, Y and Z by 

X = - 2F, Y = - 2F, Z = F[. (2-20) 

Then, from ()2-ll() . all basic gauge- invariant variables can be expressed in terms of 
these three variables as 



F = (2-21b) 



X + Y 

Inserting these expressions into (|2-12() . we obtain 

{2r/k)Et = -dtX - drZ = 0, (2-22a) 
f f 1 



(2r/k)Er = j-X - drY - ^Y + -j^dtZ. (2-22b) 



Unlike Ea, the expressions of in terms of X, Y and Z, which are given 
in Appendix are rather complicated. In order to put them into simpler forms, 
we utilize the Fourier transformation with respect to the time coordinate t. Then, 
writing the Fourier component proportional to e~"^* as 

X ^ X, Y ^Y, Z ^ wZ, (2-23) 

we find that the Einstein equations are reduced to the following equations for u; / 0: 



X' 



2/ 
Z' = X, 



r \jf 



^•y,2 



y = h^x-Y)^^z 



(2-24a) 

(2-24b) 
(2-24c) 



M 



up-r^ + K\r^ + ( n(n - \)K - n{n + l)Ar^ 



n 



DM 



>.n— 1 



+ 



2 2 

LO r 



k^f + nK"^ - (n - l)K\r'^ 



AnKM ^{n+ ifM'^ 



noj^r"^ + ( Ar 



2 {n - l)M 



r,ri— 1 



r 

Z = 0. 



n~l 



,2{n-l) 



Y 



(2-24d) 



Here, the last equation corresponds to —2r'^fE^ = 0. 

Now, we comment on the / = 1 mode. As explained in Appendix^ the above 
basic equations also hold for this mode, if we regard Ft = as a gauge condition. 
However, the number of the residual gauge degrees of freedom that remain after 
this partial gauge fixing is the same as that of the degrees of freedom of the general 
solution to H2-24() . This implies that there exist no physical degrees of freedom for 
1 = 1. For this reason, we assume that I > 2 from this point. 
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§3. Master equation 

Because the constraint ()2-24d|) is hnear, we can always reduce the system of first- 
order ODEs (|2-24() to a second-order ODE for any linear combination of X,Y and Z. 
However, the corresponding 2nd-order ODE in general has coefficients that depend 
on the frequency co in an intricate way and therefore it is not useful. In particular, 
in order to allow for its use in the stability analysis of black holes, it is desirable 
that the master equation has the form A<P = uP'<P, where A is a self-adjoint second- 
order ordinary differential operator independent of u, as does the Zerilli equation 
for the four-dimensional Schwarzschild black hole. Starting from the expression of 
the master variable for the Zerilli equation in terms of our gauge-invariant variables, 
after some trial and error, we have found that the following form of is the best: 



nZ-r{X + Y) 
Here, H \s a, function of r defined by 



H{r) = m + —n{n + (3'2) 

m = - riK, (3-3) 

2M , , 

(3-4) 



^n—l ' 



In terms of ^, (|2-24|) is reduced to 



-4Jf'4)+Vs'^ = -'^- (3-5) 



dr \ dr J 

Here, the effective potential Vs is given by 

Vsir) = ^(-)^(-) (3.6) 

with 

Q(r) = - [n^{n + 2)(n + ifx'^ - I2r?{n + l)(n - 2)mx + 4(n - 2)(n - 4)m^] y 
+n^(n + ifx^ + n(n + 1) [4(2n^ - 3n + 4)m + n(n - 2)(n - 4)(n + l)K] 
-12n [(n - 4)m + n(n + l)(n - 2)K] mx + 16m^ + AKn{n + 2)m^, (3-7) 

where 

y = \r\ (3-8) 

It is easy to see that this equation is identical to the Zerilli equatiorP for n = 2, 
K = 1 and A = and identical to the equation for even modes derived by Cardoso 
and 

LemoPJ-El forn = 2, = 0, 1 and A < 0. 

The original fundamental variables X, Y and Z are expressed in terms of the 
master variable ^ as 



f 16H^ J AH 



(3-9a) 
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Y 



Pz 

—<P - frdr^ 
4H 



(3-9b) 
(3-9c) 



Here, the coefficients Px, Qx, Py, Qy, and Pz are functions of r expressed in terms 
of X = 2M/r"-~^ and y = Ar^ as 

Pxi^) = 4(n — l)[n^(n + l)x — 2{n — 2)m\my 

+n^{n + 1)^3;^ + 2n(n + l)[2{r? +n + 2)m — n{n — 2)(n + l)-ftr]x 
-4n[(n - ll)m + n{n + l)(n - 3)K]mx + 16m^ + SKm^n^, (3-lOa) 
4(n — l)my + n{n + l)^a;^ + 2[(3n — l)m — n(n + l)ir]x 

-4Knm, (3- 10b) 



Qx{r) 
PY{r) 



4n^(n + l)(n — 3)mx + (— 8n + 16)m^]y 



[2n^(n + 1 

+n^(n - l)(n + l)^x^ + 2n(n^ - l)[4m - n(n - 2)(n + 

+An{n — l)[2>m + n{n + l)K]mx, (3- 10c) 

2[n^(n + l)x + 2m]y + n(n — l)(n + l)x^ 

-2(n- l)[m + n(n+ l)ii:]x, (3-lOd) 
[— n^(n + l)x + 2(n — 2)m\y + n(n + l)x^ 

+ [2(2n - l)m + n(n + l)(n - 2)ir]x - 2Knm. (3-lOe) 



Although we have utihzed the Fourier transformation with respect to the time 
coordinate to derive the master equation H3-5() . we can apply the inverse Fourier 
transform to it and thereby obtain the master equation in the form of a wave equation 
for a master field ^(t, r) in the two-dimensional orbit space with the coordinates (t, r). 
This master wave equation is obtained simply replacing lo by idt in (|3-5|) . This yields 



Vs_ 
f 



^ = 0, 



(3-11) 



where □ is the d'Alembertian operator in the two-dimensional orbit space with the 
metric Qabdy'^dy^. Furthermore, because the right-hand sides of ()3-9() are polynomials 
in LO, these expressions for the original gauge- invariant variables in terms of the 
master variable can be transformed into expressions for X{t,r),Y{t,r) and Z{t,r) 
in terms of <P{t, r) through the same replacement: 



X 



Y 



Z 



„n/2-l 



16H^ AH ' 



P^ 
AH 



dt$ + frdrdt<P 



If we introduce the variable fi defined by 

Q = r'^/'^H^. 



(3-12a) 
(3-12b) 
(3-12c) 

(3-13) 
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these expressions can be put into the covariant form*-' 

1 



Anr'^H . 



2H -nrf')Q + 2nrDr-DQ , (3-14a) 



Fab + {n- 2)Fgab = ^ (oaDbQ - ]pf2ga^ . (3- 14b) 

We found by symbolic computations that when these expressions are inserted 
into Ea, El, Et and Eat, the latter take forms consisting of hnear combinations of 
the master wave equation (|3-lip and its derivatives. This suggests that the set of 
equations (|3-11|) and (|3-12j) is equivalent to the Einstein equations even if the Fourier 
transforms of the gauge- invariant variables with respect to the time coordinate do not 
exist. In fact, we can confirm this by the following argument. First, by inspecting the 
procedure leading to the master equation, we find that if we take dtX(t, r), dtY(t, r) 
and Z{t, r) as basic variables and define r) by 

,^^^^^_ nZit,r)^B,Xi^rnmi^,r) ^ (3.15) 

we can obtain the time derivative of the master wave equation, 

a* (^D^ - y5>^ = 0, (3-16) 

from an algebraic combination of Ea, El, E^ and their derivatives. The expressions 
for dtX, dtY and Z in terms of dt^ corresponding to (|3-12|) can be also obtained 
with the same method. This implies that if the perturbed Einstein equations had 
a solution that could not be expressed in terms of the master wave equation ()3-ll|) 
with H3-12p . it must be static. This is consistent with the fact that original definition 
of ^ in terms of the Fourier transform, ()3-l() . becomes singular for w = 0. However, 
as we show in the next section, static solutions are also solutions to (|3-11|) with 
H3-12(l . Therefore, these equations are equivalent to the original perturbed Einstein 
equations. 

Finally, we note that in the special case M = 0, the master variable considered 
in the present paper is related to the master variable ^{s) ™ R-^f. 29) and Q in 
KIS2000 for scalar perturbations in a constant-curvature background spacetime as 

<J(^) = Q = r"/2^, (3-17) 

as is shown in Appendix O 



§4. Static perturbations 



As noted in the previous section, the master variable ^ defined by (|3-1() is ill- 
defined for (J = 0, and the derivation of the master equation for a scalar perturbation 

*■* We would like to thank an anonymous referee for recommending us to look for these covariant 
expressions. 
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given in the previous section does not apply to static perturbations. In this section, 
we show that in spite of this, the master wave equation ()3-ll() also describes static 
scalar perturbations. 

For static perturbations, from (|2-22a() . (|2-22a|) . (|A-2a|) and (|A-2b|) . we find that 
the equations Ea = 0, EJ^ = and E^'. = can be written 



Et: Z' = 0, (4- la 

L 
2/ 



f 

Er ■■¥' + ^{Y - X) = 0, (4-lb) 



El:^Z = 0, (4-lc) 



1 \ r-^ n 2r n 

IL + !11)y'- (^k + ^ + ~ 

2 r J y n 

P"-^'^' + !^/"-ar = o. {4.id) 



2r n 

Hence, Z vanishes and the basic equations are given by the following set of first-order 
ODEs for X and y: 



2(n-l) ,, , 4(n + l)A , 2f {ff^V" 

— 2 — yj ~ ^ ^ — 

n r 2j n 

2{f - K) 4(n^ - 1) , 2(n - 1) (/O^ 

2 ' '^^ J WT' 

n r 2j 

n r''' J 

Y' = ^{X-Y), (4-2b) 

Z = 0. (4-2c) 

Because this set consists of first-order ODEs for the two variables X and Y, it is 
always possible to reduce it to a single second-order ODE for any linear combination 
of X and Y. Hence, if we adopt a combination that is consistent with (|,'-{-12)) in the 
static case, it is expected that we can obtain a master equation that coincides with 
(|3-llj) without the time derivative term. In fact, there is a unique such choice, given 
by 

QY{r)X{r)-Qx{r)Y{r) 
^ ^ • 2A;2W2-i/'(r)i/(r) ' ^ ^ 

For this choice, if we express X and Y in terms of and with the help of (|4-2|) . 
we find that they are represented by H3-12a() and p-12b|) with dt^ = 0. Further, 
insertion of these expressions into (|4-2|) gives (|3-1H) without the time derivative 
term. Thus, it is found that the same master equation holds for both non-static and 
static perturbations. 
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Although it is preferable for the investigation of general cases (such as the stabil- 
ity analysis of black holes) that every perturbation be described by a single master 
equation, the potential in the master equation (|3-11() is rather complicated and is not 
always useful for analysis of static perturbations. Because the behaviour of static 
perturbations is important with regard to the issue of black hole uniqueness, it would 
be useful if we could obtain a master equation with a simpler potential. 

Now, we show that such a master equation is indeed obtained if we adopt Y as 
the master variable. First, (|4-2() leads to the following second-order ODE for Y: 



where 



Y" + aY' + I3Y = 0, (4-4) 

r'^ff'a= {n - l)(n + 2)Kx + InKXr"^ + (n - l)(n - A)x'^ 

-{r? + lln - 10)xAr2 - 2(n - 4)AV'^, (4-5a) 

r^fP = -k^ + nK + {n-2){K - f). (4-5b) 

In order to rewrite this equation in a formally self-adjoint form, we introduce the 
new variable Y defined by 

Y := SY; (4-6) 

S = /-'^e.p(/|<i,-)=-^. (4.7) 

Then, the above equation for Y is transformed into 

-4(f'4-]+yy-». (48) 



dr \ dr 



where 



Q{r)= -{2n-l){n-lfx'^ 

+2(n - l)[n(n - l)K + {br? + 23n + 6)y]x^ 

+ [n{n -2){n- ifK^ - 2{n - l)(n + 3){n'^ + lOn + 6)Ky 

+ + 28n^ + ein^ - 66n - 28)y^]x^ 
+ [4(n - l)(3n2 + lOn + 4)K'^ - 8(3n^ + Tn^ - 5n - 6)Ky 

+A{3n^ + 7n^ - 12n - 4)y^]yx 
+4n(n - 2)K^y^ - 8{n^ - 2n - 2)Ky^ + 4n(n - 2)/. (4-10) 

The dependence on k"^ of this effective potential is simpler than that of V . Further- 
more, although it is still rather complicated for A 7^ 0, it becomes quite simple for 
A = 0: 

~ 4/c2(l - x) +n(n - 2) + 2n3; - (2n - l)x2 

y = ^ — - — ^ — ^ —■ (4-11) 



13 



§5. Vector and tensor perturbations 

It was shown in KIS2000 that vector and tensor perturbations in the background 
spacetimes considered in the present paper are described by master equations of 
the wave equation type in the two-dimensional orbit space. In this section, for 
completeness, we show that these equations can also be written in the form of (|3-lljl . 

5.1. Tensor perturbations 

Perturbations of the tensor type can be expanded in terms of the tensor-type 
harmonic tensors Tij satisfying 

(A„ + 4)Tij=0, (5- la) 

= 0, DjVi = 0. (5- lb) 

The eigenvalues are all positive. They form a continuous set for K < and a 
discrete set, 

A;2 =/(/-Kn-l)-2, / = 1,2, (5-2) 

for K = 1. 

Each harmonic component of the metric perturbation is expressed by (|2-4I) . with 
fab = fa = Hl = and Sjj replaced by Tij. Ht is gauge- invariant by itself, and in 
vacuum, it satisfies the following wave equation in the two-dimensional orbit space: 



-Dr ■ DHt - 



uHt + ^Dr • DHt - ' ^ '2"" Ht = 0. (5-3) 



If we introduce the master variable ^ by 

^ = r^'/^HT, (5-4) 
this equation can be rewritten in the same form as H3-ll() : 



Vt 

□^__±^ = 0, (5-5) 



with the effective potential 



, 9 - 2) n(n 2) ^ 5 n^M 

kl + 2K + '-K ^ '-\r^ + 

^4 4 2r"-i 



(5-6) 



This potential is identical to that derived by Gibbons and HartnolPl for the more 
general case in which da"^ in (|2-14() is given by an Einstein metric, if we use the 
relation \l = k'^ + 2nK that is valid for maximally symmetric black holes, where 
Ai, is the eigenvalue of the Lichnerowicz operator. 



*' The authors thank Sean HartnoU for discussion of this poinPl 
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5.2. Vector perturbations 

Perturbations of the vector type can be expanded in terms of the vector-type 
harmonic tensors Vi satisfying 



(A„ + 4)Vi = 0, 
DjY^ = 0. 



(5-7a) 
(5.7b) 



As in the case of the tensor-type harmonics, the eigenvalues ky are all positive, 
forming a continuous set for K < and a discrete set for K = 1. However, this 
discrete spectrum is shifted by unity with respect to that in the tensor case: 



fey 



/(/+n- 1) - 1, l = l,2,-- - . 



(5-8) 



The expression for the metric perturbation is now given by H2-4() with Sj replaced 
by Vj and Sij by 

'Vij = -:r]-{DiYj+DjYi), (5-9) 

2fcy 

and the non- vanishing harmonic coefficients are fa and Hx- For /c^ > (n — 1)K, 



Fa = fa + ^DaHr 
ky 



(5-10) 



is gauge invariant and represents a natural basic variable. In our vacuum background, 
F"" is expressed in terms of a field i7(t, r) in the two-dimensional orbit space satisfying 
the wave equation 



n 

UQ - -Dr ■ DQ 

r 



as 



„n — 1 r?a 



where Cab is the Levi-Civita tensor of the two-dimensional orbit space. 
If we introduce the master variable ^ by 



(|5-11|) takes the form 



with 



Vvir) 



kl + K + 



Vv 

n<p- -^4> = 0, 



n{n-2)K n{n - 2) 



\r^ 



3n^M 



(5-11) 
(5-12) 

(5-13) 
(5-14) 

(5-15) 



This equation is identical to the Regge- Wheeler equation for n = 2, K = 1 and 
A = and identical to the equation for odd modes derived by Cardoso and Lemos 
forn = 2, K = 0, 1 and A < 
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For ky = {n—l)K, i.e., K = 1 and / = 1, Yij vanishes, and the perturbation vari- 
able Ht loses meaning. In this case, the following becomes a basic gauge-invariant 
variable: 

F;;) = rDaih/r) - rDbifa/r). (5-16) 

In our vacuum background, from the Einstein equations, it follows that is ex- 
pressed as 



(1) _ Le^ 



where L is an arbitrary constant. This simply represents a rotational perturbation of 
the black hole, and in the spherically symmetric case, L corresponds to the angular- 
momentum parameter in the Myers- Perry solution.^'" 



§6. Discussion 



In this paper, we have derived master equations of the wave equation type that 
describe gravitational perturbations of maximally symmetric black hole spacetimes in 
higher dimensions. For scalar-type and vector-type perturbations, respectively, they 
represent extensions of the Zerilli equation for polar perturbations and of the Regge- 
Wheeler equation for axial perturbations of the four-dimensional Schwarzschild black 
hole to higher dimensions as well as to the case of a non- vanishing cosmological con- 
stant and to quasi-black hole spacetimes in which constant-time sections of equipo- 
tential surfaces have non-positive sectional curvatures. Hence, these master equa- 
tions are expected to be useful in a wide variety of higher-dimensional gravity prob- 
lems. 

Our formulation also gives extensions of the Zerilli and the Regge- Wheeler for- 
malisms to the four-dimensional Schwarzschild-de Sitter and Schwarzschild-anti- 
de Sitter backgrounds, which coincide with those given by Cardoso and Lemos 
recently.^ Here, we show that we can prove the stability of four-dimensional 
Schwarzschild-de Sitter and Schwarzschild-anti-de Sitter black holes using our for- 
mulation. For a spherically symmetric spacetime of four dimensions, i.e., for n = 2 
and K = 1, Q in (|3-7|1 becomes 

Q = 288xV(^) + 432x^ + (288 + lUlj + 720/2)x^ 

+(^2 + 1)^(^2 + 4)2[482; + 16(^2 + 3)(/2 + 2)], (6-1) 

where I2 = I — 2. Hence, the effective potential Vs in (|3-5|) is positive definite. This 
suggests that in (|3-5|) is always positive and that the black hole is stable with 
respect to scalar perturbations. In fact, this is actually the case for a non-extremal 
Schwarzschild-de Sitter black hole, if we consider the region bounded by the black 
hole horizon and the cosmological horizon, because the range of r* = J dr/f is 
(—00, +00), and becomes convex for lj^ < (changing the sign of ^ if necessary). 

This simple argument does not hold for the Schwarzschild-anti-de Sitter black 
hole, as r^, has a finite limit for r ^ 00 in this case. This is closely related to the fact 
that the differential operator on the left-hand side of ()3-5|) does not have a unique 
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self-adjoint extension in the L^-space with respect to the inner product 

{<Pi,<P2)l^ = J dn^lir)^2ir). (6-2) 

In this case, we must impose a boundary condition at r = oo to make the problem 
well-posed. The most natural choice is to require that ^ vanish at r = oo, if we are 
concerned with the local stability of the black hole. Under this boundary condition, 
we can apply the same reasoning as for the Schwarzschild-de Sitter black hole, and 
we thereby find that the Schwarzschild-anti-de Sitter black hole is stable. 

We can also undertake a similar analysis for vector perturbations of these four- 
dimensional black holes, although the situation in this case is slightly subtle. Here, 
the effective potential becomes 

Vy = l(-— + 6 + h{h + 5)") . (6-3) 

For yl > 0, this is positive definite for / > 2 outside the horizon from 1 — 2M/r = 
/ + Ar^ > 0. Hence, the argument given above again applies. In contrast, for /I < 0, 
Vy becomes negative near the horizon for large |A|, and the above argument does 
not apply. Nevertheless, we can prove using the somewhat sophisticated argument 
given in Ref. 41) that there exists no unstable mode in this case as well. 

From the above considerations, we have established the perturbative stability 
of spherically symmetric black holes in four dimensions, irrespective of the sign 
of the cosmological constant. The above argument also proves the uniqueness of 
these black holes in a perturbative sense, because there should exist a regular static 
solution to the perturbation equation with Z > 1 if there exists a continuous family of 
regular static black hole solutions with the same mass that contains the spherically 
symmetric solutions. 

In spacetimes of dimension greater than four, the stability of spherically sym- 
metric black holes has not yet been investigated, even in the asymptotically flat case. 
It seems that our formulation would be useful in the investigation of this problem. 
For example, using our formulation, we can prove the stability of higher-dimensional 
Schwarzschild black holes, as will be shown in a separate paper.SS This is not as 
trivial as in the four-dimensional case discussed above, because the effective poten- 
tials in the master equations are not positive definite in general for vector and scalar 
perturbations. Nevertheless, it can be shown that the differential operators on the 
left-hand sides of H3-5() and (|5-14() have unique positive self-adjoint extensions. A 
similar technique might be used to analyse the asymptotically de Sitter and anti-de 
Sitter cases in higher dimensions as well, although the behaviour of the potential in 
these cases is much more complicated. This problem is now under investigation. 

Next, we comment on the relation between the master variable for a scalar 
perturbation and that for a vector perturbation, which we here denote <Ps and ?^y, 
respectively. In the four-dimensional Schwarzschild background with yl = and 
K = 1, Chandrasekhar and Detweiler found that they are related by the simple 
relation 

= P^v + q^'v, (6-4) 
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where p and q are appropriate functions of r that are independent of the frequency 
w if and only if the corresponding effective potentials are expressed in terms of a 
single function F of r independent of uj as 

dF 

Vs. Vv = ±f^ + F^ + cF, (6-5) 

with some constant c. This relation with c = k'^{k'^ — 2K) / {3M) holds for any 
values of A and K in four dimensions, as pointed out by Cardoso and LemosP^'^ 
However, we have found that there does not exist a function F satisfying ()6-5() . 
and hence, there is no simple relation of the form (|6-4j) between the two types of 
perturbations in higher dimensions. This result may be closely related to the fact 
that a tensor-type perturbation appears as a new mode in higher dimensions, and it 
implies that the spectral analysis of gravitational perturbations may not reduce to 
that for a vector-type perturbation with a simpler effective potential, unlike in four 
dimensions. 

Finally, we comment on an extension of our formulation. As mentioned in §2, 
the background metric (|2-14() with (|2-15() satisfies the Einstein equations even if the 
metric da'^ for /C" is replaced by an arbitrary Einstein metric satisfying 

Rij = in-l)Kjij. (6-6) 

In this case, the Weyl curvature C^jki of /C" no longer vanishes, and it provides a non- 
trivial background tensor distinct from 'jij. However, because the Weyl curvature is 
trace-free and of second order with respect to the spatial derivative, it can couple 
only to a tensor perturbation in the linear theory. Hence, the same perturbation 
equations hold in this generalized case for scalar and vector perturbations. As a 
consequence, the master equations for these types of perturbations derived in this 
paper are also valid in the case that /C" is Einstein, if we replace the eigenvalues kg 
and ky by the corresponding ones for the Laplacian in the Einstein space. Further, 
for tensor perturbations, C^jki appears in the perturbation equation only through 
the combination called the Lichnerowicz operator, 

{ALh)ij = -Anhij - 2Rikjih^^ + 2(n - l)Khij, (6-7) 

where hjj is a perturbation of the metric 5gij of /C", as shown by Gibbons and 
Hartnoll.!^ Hence, if we expand the tensor perturbation in terms of the eigentensors 
with respect to A/, instead of A„ and replace the eigenvalue A/, of by A;|n + 2nK. 
we obtain the same equation for tensor perturbations, as mentioned in §5.1. More 
detailed explanation of this extension will be given in a separate paper. 
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Appendix A 

The expressions of Eat, Ea, and Et 



We have the following: 



2E. 



ab 



-UFab + DaD,F§ + DhD.F^ + n {-D.Fab + DaF,b + DbF,a) 



+ 










n J 


-2n DaDbF^ 








r 






— D r 




r 



2n 



n(n — 1] 



D'rD'^r ] F^a 



2n(n +1) , ^k'^ -nK 
-2naF ^ ^-Dr ■ DF + 2 n - 1) ^ F 



9ab, 



2El 



^ Dbir'^-^F^) + rDa (-F^ ) + 2(n - 1)Z),F, 



r-n— 2 



(A-la) 
(A-lb) 



-DaDbF"^ - — ^^D^rD^Fab - '-^-^ i {n - 2)D''rD\ + 2rD''D\ ] Fab 



n — 1 



n 



2n(n-l)^ ^„ 2(n - l)(n - 2)(k'^ - nK) ^ ,^ , 
+2(n-l)nFH ^ '-Dr-DF--^ '-^ 2^ ^F, (A-lc) 



2r\ 



-2{n-2)F-Fl 



(A-ld) 



Here, A is the cosmological constant and i?*^^) is the scalar curvature of the two- 
dimensional metric Qab- 

For the metric (|2-14() . the quantities E^ are expressed in terms of the variables 
X, y and Z defined in TTM as 



2Fr 



^-r-^-2(n + l)A 

+fdtdrY +{^-^-\dtY 



z 



+fdtdrX 



(A-2a) 



2f: 



1 f 

d?x - ^—drX 



f 
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n 



2r 



n 



X 



1 



f 



2 r 



+ 



'K-f 2(n2-l) 3n-2 , n - 1 „ - nK' 
— / / 19 — 



n 



2r 



n 



2n^ 
rf 



Y 



(A-2b) 



2Ei = -fdlX + {^f - l^^j drX 

(2n - 3)/ 2{n'^ - 1) n - 2 n - 1 k 



-K 



+ 



A + 



n 



2r 



n-1 n-3, 2(n + l) ^ , {n-2)f' f" 

o 

n 2r n 



n 



Y. 
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-/' + — r - ^ 



X 



(A-2c) 



Here, A 
El: 



2A 



n{n+l) 



. For completeness, we also give the corresponding expression for 



El = ^d^X + ^drX 



2/ 

+ 



4 



(n-l){n-2){f -K) ^ (6n - 4 - n^)f ^ T 
2nr2 4nr 2n 

3/' , / 



+ 



{n-l){n-2){f - K) (-n^ + 2n - 4)/' 



2nr'^ 



+ 



4nr 



2n 



1 /' 



r/ 2/^ 



1 



+ {—r-:hj]dtZ + -dtdrZ. 



f 



(A-3) 



Appendix B 

The I = 1 mode in the spherically symmetric case 



For the mode with Z = 1 in the spherically symmetric case, the perturbation 
variable Ht and, correspondingly, the component Et of the Einstein equations do 
not exist, since Sij = 0. However, by introducing the variables Fab and F using (|2-7|) 
with Ht = 0, we can recover the equation Ej- = as a gauge condition. To see 
this, first note that Fab and F are no longer gauge invariant, and from the equations 
(48)~(52) in KIS2000, they transform under the gauge transformation 



Sy" = T^s, 5z' = m\ 



(B-l) 
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as 



SFab = {rDaDbL + DarDbL + DbrDaL) 



6F 



--{Dr-DL + ^ ^—L 

k \ nr 



In particular, Et transforms as 

1 



(B-2a) 
(B-2b) 

Prom this, it follows that we can always make Et vanish by choosing an appropriate 
gauge. This requirement does not fix the gauge completely, and there remain resid- 
ual gauge degrees of freedom parametrized by two arbitrary functions of the space 
coordinate r in the two-dimensional orbit spacetime 

Appendix C 

The relation to the master variable in the case M = 



5Et = ^D, (r-D'^L) + 



The master variable ^(5) in Ref. ISHJ and Q (= ^(s)) in KIS2000 for scalar 
perturbations are related to the gauge-invariant variables by 



In 



Fab = DaDbf2 - 

and they satisfy the wave equation 

71 

UQ - -Dr ■ DQ - 

r 



n I V n 



k^ — uK , X , 1 ^ 
2 + (n-2)A) f? = 0. 



Z is expressed in terms of fi as 

Z = D'DtQ = f/^dr 
Also, for M = 0, H3-12c|) becomes 

Z = f/'dr 

Hence, we have 

Because it is easily checked that f2 = r"/^'^ satisfies ()(]-2|) . this relation implies that 
C{t)y/f also satisfies this equation. This condition is expressed as 



(C-la) 
(C-lb) 

(C-2) 

(C-3) 

(C-4) 
(C-5) 



C 



{m + K)X 



mK 



C. 



(C-6) 
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From this, it follows that K = if C ^ 0. However, for K = 0, the expression for X 
in terms of f2, 

X = -jDtDtn-ia + \)n, (c-7) 



is identical to (|3-12ap only when 



C = --n{n - 1)AV2c. 



(C-8) 



Because A 7^ for = 0, this implies that C = 0. Hence, 17 is related to ^ by 
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